Magnetic Properties for the One-Dimensional Multicomponent Spin-Gap System 



Akira Kawaguchi, Tatsuya Fujii and Norio Kawakami 
Department of Applied Physics, Osaka University, Suita, Osaka 565-0871, Japan 

(February 1, 2008) 



O 
O 

(N 

a 

ON 



i 

CO 

c3 



o 
o 



> 
m 

^1" 



c3 



i 

c 

o 
o 



X 
S3 



Magnetic properties for the one- dimensional multicompo- 
nent quantum spin system with the excitation gap are studied 
based on the integrable spin model introduced by Bariev et 
al. By exactly computing the magnetization, we show how 
the characteristic structure with plateaus and cusps appears 
in the magnetization process. To study low-energy dynam- 
ics of the system, we apply the finite-size scaling analysis 
to the excitation spectrum, and thereby evaluate the power- 
law exponent as well as the enhancement factor for the low- 
temperature NMR relaxation rate 1/Ti. We discuss the criti- 
cal properties of 1/Ti around plateaus and cusps in the mag- 
netization curve. 

PACS numbers: 75.10.Jm, 75.40. Gb 



I. INTRODUCTION 

Quantum phase transitions of one-dimensional (ID) 
spin systems in a magnetic field have attracted much in- 
terest recently. One of the most remarkable phenomena, 
is the plateau formation in the magnetization process,ETu 
which is caused by a field-induced spin gap. Also, the 
closure of the spin gap in the Haldane system and the 
ladder system by the magnetic field has been providing 
an interesting subjectjjij which shares common impor- 
tant physics with the plateau formation. 

Another hot topic in ID spin systems has been con- 
cerned with the multicomponent version of the spin 
model or the spin-orbital model, since it was recog- 
nized that the orbital degrees of freedom give rise to a 
variety of interesting phenomena. Experimentally, ID 
correlated electron systems-, at quarter-filling, such as 
Na 2 Ti 2 Sb 2 0el and NaV 2 5 H have been studied exten- 
sively, for whict,tie orbital degrees of freedom play an 
important role.t3~t3ln this context, the ID SUI^Wiass- 
less spin mpdeL has been studied numericalljOEfl and 
analyticallyLaiij as the simplest spin-orbital model, and 
this type of analysis has been extended to .a, perturbed 
SU(4) model to study the gap formation.E3cil 

According to the above studies, it may be interesting 
to investigate how the orbital degeneracy affects field- 
induced quantum phase transitions such as the plateau 
formation in ID spin systems. To address this problem, 
in this paper, we investigate static and dynamical prop- 
erties for the ID multicomponent spin-gap system in a 
magnetic field. For this purpose, we employ the ID mul- 
ticomponent anisotropic Heisenberg model, which was 
exactly solved by Bariev et al.rQ In particular, we focus 
our attention on the critical regions where the spin-gap 



disappears, or the magnetization plateau (or cusp) is gen- 
erated in the presence of a magnetic field. To study the 
low-energy dynamics of the system, we further discuss 
the NMR relaxation rate 1/Ti around the field-iaduced 
quantum phase transitions at low temperatures. BjEj 

This paper is organized as follows. In Sect. II, we 
exactly calculate the magnetisation curve for the ID 
anisotropic Heisenberg model^a and determine the mag- 
netic phase diagram. We clarify how the plateaus and 
cusps are generated in the magnetization curve— In 
Sect. Ill, by exploiting finite-size scaling techniquesLJ we 
discuss the low-temperature behavior of the NMR relax- 
ation rate 1/Ti with the emphasis on the critical proper- 
ties around plateaus and cusps. It is shown that the en- 
hancement factor plays a significant role for 1/Ti around 
such critical regions. A brief summary is given in Sect. IV. 



II. PLATEAUS AND CUSPS IN THE 
MAGNETIZATION 

In this section, we study the magnetization process 
for the anisotropic Heisenberg model with a two-band 
structure in ID, which is defined-by the deformed version 
of the SU(4) Heisenberg model f5 

Hq = JY^ {4m c irn'4 +lm ,C i+lm 
i m^m' 

-exp[sign(m' - m)$]n im n i+ i m >}, (1) 

where it is prohibited that more than one electron oc- 
cupies each site. Here, c\ m is the creation operator for 
electrons with spin (f , |) and orbital (1, 2) degrees of free- 
dom. We shall use the indices m = 1, 2, 3, 4 for these four 
states so that the corresponding number of electrons N m 
should satisfy, Ni > -/V 2 > A3 > A4 in a magnetic field. 
In the isotropic case ($ = 0), the Hamiltonian_/|l|) is re- 
duced to the SU(4) massless Heisenberg model,E3 whereas 
for any other cases it is massive for all excitations at zero 
magnetic fieldE3 Although the above integrable model is 
rather special in its appearance, we think that its char- 
acteristic behavior in a magnetic field, such as the forma- 
tion of plateaus and cusps, should capture some generic 
properties expected for multicomponent spin-gap models 
in ID. 

The exact solution of the model ([!]) with periodic 
boundary conditions was obtained by Bariev et al.cB The 
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mations obtained for rapidities (a = 1, 2, 3) 
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toclf - N®x{Xj)5 a i = E ©i(^ - \ Q+S ) 

s=±l i=l 

-^9 2 (A^-Af), (2) 



A. Spin-orbit multiplet case 

We first assume that the Zeeman splitting is given by 
g^sMzH with M z — ±1/2, ±7/2, where 7 ranges from 
to 1 . The corresponding Hamiltonian is 



with M a = E m =a N m (M = Mi = 0), where 9 n (A) = 
2 tan -1 (cosh tan -|) . The quantum number Ij^ which 
specifies elementary excitations is subject to the con- 
straints, 
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(Mi + M 2 ) mod f, 
--(M1-M2 + M3) mod f, 
'-{Mi - M 3 ) mod f. 



(3) 



In the thermodynamic limit at zero temperature, the al- 
gebraic equations (|^) are converted to the linear integral 
equations for the distribution functions of rapidities a a 
and for the dressed energies e a , respectively ,t3o 

3 r+A °dA; 

K ai (\ a — A 7 )e 7 (A 7 ), 



£ Q (A G 



E 

7=1 



2?T 



(4) 



^(Aq.) = er Q (A Q ) 



3 

E 

7=1 



'^K (A 

Z7T 



where K a 



A'. 1 



A> 7 (A 7 ), (5) 



a, 7+1 



J^n(A) = sinh(n$)/[cosh(n$) -cosA]. 
distribution function is given by 

^a(Al) = ^"^(Ai)^!, 
Z7T 



), and 
Here, the bare 



^a+1,7 



where n i a denotes the particle number for electrons 
with M z = a. Note that the case of 7 = 1/3 corre- 
sponds to the Zeeman splitting for the spin-orbit multi- 
plet J z = ±3/2, ±1/2 (M z = \ J Z ), while for 7 = 1, the 
Zeeman effect acts only on the spin sector, so that the 
four-fold multiplet splits into two orbitally-degenerate 
levels in magnetic fields. The Zeeman term (|?]) may be 
applied to more generic cases including these two impor- 
tant cases. For example, in some Ce compounds such 
as CeB 6 the four-fold ground multiplet (r 8 multiplet) is 
generated in a cubic crystalline field under strong spin- 
orbit coupling. In such cases, the value of 7 is known 
to depend on the strength of the crystalline filed. We 
thus regard 7 as a continuously changing parameter in 
the following discussions. We define the magnetization 
in terms of the distribution functions, 

m (A) = {(JVi/a ~ N-1/2) +7(^7/2 " #-7/2)} l L 

= l- 2 7 / cr 2 (A 2 )dA 2 

J-X° 

( + A° +A° 1 

-(1 - 7) I J J tri(Ai)dAi + J J a 3 (A 3 )dA 3 I , (8) 



(6) where we have set the number of electrons as Ni 



N 1/2 ,N 2 



-sinh$X 1 (Ai) + i(l- 7 )/i 
7/1 



and e a will be explicitly given separately for each case. tively. Here e° in (||) is given as 
The cut-off parameters ("Fermi level" for elementary ex- 
citations) A" (0 < A" < tt) are determined by minimizing 
the free energy-ibr a given magnetic field, resulting in 
the condition,c3a e Q (±A° ) = 0. For zero magnetic field, 
they are reduced to A° = X 2 = A3 = tt. Since the dressed 
energies satisfy £i(±Aj ) = tt) = £ 2 (±A2 = tt) = e 3 (±A 3 = 
tt) < 0, the system has three kinds of spin-gap excita- 
tions, which are degenerate at zero magnetic field. 

As mentioned above, one of the most remarkable phe- 
nomena in ID spirusystems is the plateau formation in 
the magnetization,Eri3 which is regarded as a quantum 
phase transition driven by the. magnetic field. Also, the 
middle-field cusp singularity^ is another typical exam- 
ple of the field-induced quantum phase transition. In the 
following, we show that the multicomponent model (|l|) 
possesses a rich phase diagram in a magnetic field, where 
the magnetization has plateaus and/or cusps depending 
on the anisotropy as well as the Zeeman splitting. To 
demonstrate these things, we deal with two typical cases 
for which the energy splitting due to the magnetic field 
appears in different ways. 



N~ / / 2 ,N 3 = iV_ 7/2 and N 4 = N_ 1/2 , respec- 



(9) 



We have performed the calculation by iterating eq.(Q) 
to determine the cutoff parameter A 7 for a given mag- 
netic field, and then have evaluated the magnetization 
by solving eq.(||). Let us start with the spin-orbit level 
scheme for the multiplet specified by the total angular 
momentum. In Figs. 0(a) and (b), the magnetization 
curve for 7 = 1/3 is shown by choosing two different val- 
ues of the anisotropy, $ = 1.7 and <& = 1.5. When the 
magnetic field is increased from zero, the magnetization 
starts to increase at the finite critical field with a square- 
root dependence. For $ = 1.5, the m — 1/3 plateau 
appears, which is followed by two cusps in higher mag- 
netic fields. On the other hand for $ = 1.7, two plateaus 
appear at m = 1/3 and m — 2/3, and the cusp struc- 
tures are not generated. We note that the formation 
of plateaus in the present model is different from the 
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FIG. 1. The magnetization curve for (a) $ = 1.7 and 
(b)$ = 1.5: We denote the quantum phase transition points 
by the arrows. In particular, the positions of cusps are shown 
by the arrows with broken line. The saturated magnetization 
is m — 1. 



standard fractional quantization discussed by Oshikawa 
et alB, for which the plateau is stabilized by the spin- 
gap state that has magnetic order commensurate with 
the underlying lattice structure. On the other hand, in 
the present case, the plateau is generated when one of 
the multicomponent spin modes becomes massive, and 
the corresponding ground state is still in a disordered 
state which is not accompanied by the commensurate 
spin arrangement. 

In Fig. the h — $ phase diagram is shown for the 
case of spin-orbit level scheme with 7 = 1/3. Note that a 
massless mode exists in the regions I, II and III, in which 
the magnetization increases continuously with the field. 
In the area where two of I, II and III overlap with each 
other, there exist two massless modes. The cusp struc- 
ture appears just at the boundary where these massless 
regions overlap. On the other hand, the regions of (a)- 
(d) outside of the massless regions denote the spin-gap 
phases, where the magnetization exhibits plateaus quan- 
tized as (a)m = 0, (b)m = 1/3, (c)m = 2/3, (d)m = 1, 
respectively. Note that -the $ = line corresponds to the 
isotropic SU(4) model Jl3 which is known to have only 
two cusps without plateaus £3 The plateaus appear at 
to = 1/3 and m = 2/3 for sufficiently large $. 
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FIG. 2. The h — $ phase diagram for the model with 
(0) for 7 = 1/3: The regions I, II and III denote the areas 
which are respectively sandwiched by two dash-dotted lines, 
two broken lines, and two solid lines, where each elementary 
excitation, ei, £2, or £3 is massless. Magnetization plateaus 
appear in the spin-gap phases outside of the above three re- 
gions: (a) m = 0, (b) m = 1/3, (c) m = 2/3, (d) m — 1. 
Magnetization cusps are observed at the boundary where two 
kinds of massless regions overlap. 



For reference, we show in Fig. || the h — 7 phase dia- 
gram for <f> = 1.5 to see how the scheme of the Zeeman 
splitting affects the magnetization curve. Recall that for 
the case of 7 = 1, the four- fold degenerate states split 
into two-fold orbitally degenerate states in a magnetic 
field, for which the plateau-structure is absent. If the 
value of the parameter 7 deviates from 7 = 1/3,1, the 
plateau appears at m = (1 + 7)/2 in the region (c). 



B. Multiplet in a crystal field 

We give another example of the magnetic phase dia- 
gram by employing a slightly different level scheme for 
which the orbital-splitting A due to the crystal field ex- 
ists at zero magnetic field, and these two distinct levels 
further split in the presence of the magnetic field h. The 
corresponding Hamiltonian is 



H 



(B) 
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~/i(n iiT - n a ) - A(n»,i 



n»,2J 



(10) 



where nj = rc-ij+n^f andrti = nif+nij. The calculation 
of the magnetization is performed straightforwardly in 
the way outlined above. For h < 2A, we set the number 
of electrons as Ni = N^, N 2 = N lt , N s = N 2 -\ and iV 4 = 
N21, respectively. In this case, e° a is given by 
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FIG. 3. The h — 7 phase diagram for the model (|l|) with 
(0) at <!> = 1.5 Each region is specified following the way in 
Figl 



On the other hand, 
electrons Ni — N^, 
for which e° reads 



for 

No : 



h > 2 A, we set the number of 



N 2h N 3 



N±\ and N 4 = N : 
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■sinh$Xi(Ai) + 2 A 
h-2A 
2A 



2,1 > 



(12) 



In Fig. ^, we show the ft, — A phase diagram derived for 
the anisotropy parameter $ = 1.5. Note that the phase 
diagram has two different regions labeled by II. This is 
because both the magnetic field and the crystal field play 
a similar role for the formation of the spin gap. The areas 
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FIG. 4. The h - A phase diagram for Eqs. (B) and mm at 
$ = 1.5. The meaning of massless regions I, II, III is same as 
in Fig.^. In the areas outside of these regions, the system is 
in a spin-gap phase. 

outside of the regions I, II and III are all massive with 
the spin gap, as mentioned above. The mechanism for the 
formation of plateaus and cusps is as discussed before. 

As summarized above, the multicomponent spin sys- 
tem with the spin gap exhibits interesting features such 



as plateaus and cusps in the magnetization curve, de- 
pending on the anisotropy parameter. 



III. NMR RELAXATION RATE 

We have so far treated the static magnetization, and 
have shown that the quantum phase transitions are in- 
duced by the magnetic field, which are accompanied by 
the plateau and cusp singularities. It is then interesting 
to observe how dynamical quantities behave around such 
field-induced quantum phase transitions. To address this 
problem, we here study the NMR relaxation rate 1 /T\ in 
magnetic fields with the emphasis on the behavior around 
plateaus and cusps. 

The NMR relaxation rate for the multicomponent 
Tomaaaga-Luttinger liquid is given at low temperatures 

bym 
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lim 



2k B T fdk 



h 2 u 



2tt 

J2B a T a (h)T"" 



A 2 (k)Imx(k,uj) 



where 



T a (h) 



n 



2.T„ 



(13) 



(14) 



and B a is a constant related to the hyperfine coupling, 
which is assumed to be independent of the temperature 
and the magnetic field. The above power-law dependence 
in temperature is typical for_correlation functions in the 
Tomonaga-Luttinger liquid.Lj Since we are now dealing 
with the multicomponent system, there are several kinds 
of relaxation processes, providing different critical expo- 
nents r] a , where the index a classifies elementary exci- 
tations, e.g. 2fc^-current excitation, etc. To determine 
•q a , we should further incorporate, the interference be- 
tween different excitation modesJ^J which are explicitly 
expressed as rj a = 2z QM , where x afi is the scaling di- 
mension originating from the interference between a and 
/i modes. This formula can be straightforwardly obtained 
by analyzing the excitation spectrum in terms o£fpjtc- 
size scaling techniques in conformal field theory.LZliiil 

In the above formula, we have explicitly written down 
the factor (|l4|) depending on the renormalized velocities 

, which have been usually neglected in the discussion of 
the temperature dependence of 1/Tj. However, in order 
to study the field-dependence of 1/Ti, we should take into 
account this factor, since it plays a crucial role around 
the critical points, as will be shown momentarily We 
will refer to this factor as the enhancement factor in the 
following discussions. 

Note that if the NMR is done on the nucleus of lattice 
spin, the relaxation occurs through a contact interaction 
and thereby 1 /T\ depends only on the transverse suscep- 
tibility X-L- O n the other hand, if it is done on other 
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FIG. 5. Critical exponents around the m = 1/3 plateau in 
Fig. ^(a). The solid lines denote the most relevant critical ex- 
ponents for the transverse spin susceptibility while the broken 
lines that for the longitudinal susceptibility. 



critical exponents are smaller (larger) than unity, 1 JT\ di- 
verges (vanishes) with power-law dependence as the tem- 
perature decreases (see (|l3|)). Similar behaviors can be 
found for the one-comppaent spin-gap system studied by 
Chitra and Giamarchi.Ea We here pay our attention to 
the enhancement factor T a at sufficiently low tempera- 
tures. In Fig. ^, we show T±_ for the most relevant sector 
in the transverse susceptibility as a function of the mag- 
netic field. We note that F\\ and T±_ exhibit essentially 
the same critical behaviors. It is seen that 1/Ti is consid- 
erably enhanced in the vicinity of the plateau according 
to the dramatic renormalization of the velocities in the 
vicinity of the spin-gap phase. The enhancement of the 
relaxation rate is observed whenever the system is in a 
massless phase close to the spin-gap phase. For example, 
the relaxation rate is enhanced when the magnetization 
decreases to zero before forming the zero- held spin gap. 



neighboring nuclei, the relaxation is through dipolar in- 
teractions, and thus 1/Ti depends not only on but 
also on the longitudinal susceptibility x\\- ln order to 
treat generic cases, we take into account the contribu- 
tions both from x± an d X\\- I n the following, we discuss 
the critical properties of 1/Ti around plateaus and cusps. 



A. Around the plateau 

We start our discussion with the massless regions close 
to the magnetization plateau. For instance, we look at 
the critical region beside the m — 1/3 plateau in Fig. 
|l|(a), where only one massless mode exists both in the 
left-hand side (III) and the right-hand side (II). In the 
region III, the scaling dimension x is obtained in a stan- 
dard form via finite-size scaling analysis £3 



x = —~ AM 2 
4£f 



£3 A-Er 



(15) 



B. Around the cusp 

Let us now observe what happens for 1/Ti around the 
cusp structure in the magnetization. We first look at 
the cusp at h ~ 2.1 in Fig. |l|(b). The critical behavior 
in the left-hand side of the cusp is described in terms 
of a massless mode. Therefore, this region can be ana- 
lyzed following the way done for the plateau case. We 
thus have the enhancement of 1/Ti near the cusp. On 
the other hand, on the right-hand side of the cusp, two 
massless modes appear which need a different treatment 
for their critical properties. According to the conformal 
field theory analysis of multicomponentpCases, the scaling 
dimension is given by the matrix form,E31iJ 



1 



x = -AM 1 (£ 

4 



(£- L )AM + AD T $, T £AD, (17) 



where the quantum numbers in the vector representation 
are AM = (AM U AM 2 ) T and AD = (AD 1 ,AD 2 ) T , 
which are subject to the selection rules, 



where so-called dressed charget 
determined by 



I £3 = hm A A o £3 (A 3 ) is 



£3 (A 3 ) = 1 - 



+ A 3 HA' 

_3tf a (A3-A^ 3 (A£). (16) 

\0 ZTT 



Here AM (AD) is the quantum number which changes 
the number of spinons (carries the 2kp current). For 
the longitudinal spin susceptibility x\\ the most relevant 
critical exponent is given by 7-711 — 2x — 2£f by choosing 
AM = and AD = 1, while for the transverse spin 
susceptibility \±, it is given by r]± = for AM = 1 

and AD = 0. The critical exponents in region II can be 
obtained via a similar analysis. 

In Fig. [5], we show the computed critical exponents as 
a function of the magnetic field h. Note that when the 
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FIG. 6. The enhancement factor T± around the plateau, 
i = 1/3, in the magnetization curve. 
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FIG. 7. The critical exponents around the cusp. As an 
example, we show them for the cusp of h -2.1 in Fig. 0(b). 



AL>i = -AM 2 mod 1, AD 2 



-AMi mod 1. (18) 



Here the 2x2 dressed charge matrixEj £ has the elements 
= Capi^) which are determined by 



r +x ° d\' 

+ E / if K ^ X ? ~ ^W^), (19) 

7=1,2 77 

where a,/3,7 = 1,2, and Kn(X) = K 22 (\) = -K 2 (\), 
K 12 (A) = K 21 (A) -Ai(A). 

According to ref.Ej, by properly choosing the set of 
quantum numbers (AMi, AM 2 , ADi, AD 2 ), we can de- 
termine the most relevant critical exponents. For Xll, 
we set (0,0,-1,0), (0,0,-1,-1) and (0,0,0,-1), be- 
cause there is no change in the z-component of spin 
(AMi = AM 2 = 0), while for x±, we set (-1,0,0,+i), 
(-1,-1, and (0,-1, -±,0). In Fig. |. the most 
relevant critical exponents obtained for x± arL d Xll are 
shown as a function of the magnetic field. It is seen that 
the critical exponent r]^ is smaller than T]±_ in a wide range 
of magnetic fields, which implies that x\\ is dominant for 
the relaxation. As mentioned before, when the dipolar 
interactions are weak, Au should be nearly zero, so that 
the divergent behavior in 1/Ti may be masked at low 
temperatures in that case. On the other hand, the di- 
vergent behavior should be always observed around the 
plateaus, at least either on the right or left region. Fig. || 
shows the results for the enhancement factor T± around 
the cusp in the magnetization. In contrast to the case of 
the plateau, the enhancement of 1/Ti appears on either 
side of the cusp, because r Q is enhanced when one of the 
massless modes becomes massive as the magnetic field 
is changed. For example, in the case of $ — in Fig. 
||, we have several quantum phase transition points at 
h ~ 0.54, 1.4, 6.0, for which we can observe the enhance- 
ment of 1/Ti in the left-hand side of these points. 

We have so far mentioned the enhancement of 1/Ti 
around plateaus and cusps in the magnetization. It is 



noted here that such enhancement of 1/Ti is induced 
not only by the magnetic field but also by the orbital 
splitting. For instance, in Fig. ||, the relaxation rate for 
h = is enhanced near the both edges of region II, when 
the orbital splitting A is continuously changed. 



IV. SUMMARY 

We have studied the magnetization process as well as 
the spin dynamics through the NMR relaxation rate for 
the ID multicomponcnt spin-gap system. By using the 
solvable model of Bariev et al, we have found that the 
plateaus and the cusps are generated in the magnetiza- 
tion curve. We have also studied the NMR relaxation 
rate 1/Ti near plateaus and cusps in the magnetization. 
For the critical region around the plateau, 1/Ti diverges 
at low temperatures at least on either side of the plateau, 
whereas for the critical regions around the cusp, it de- 
pends on the values of the hyperfine coupling constants 
whether such a divergent behavior appears or not. It has 
also been found that the relaxation rate at low tempera- 
tures is enhanced near the critical value of the magnetic 
field. This is not only the case for the multicomponent 
model but also for ordinary spin chain models in mag- 
netic fields. The enhancement appears at sufficiently low 
temperatures, so that it may be observed in rather ideal 
ID spin systems for which the long-range order due to 
3D dimensionality is suppressed even at low tempera- 
tures. Also, to extract this enhancement effect from the 
experimental data, careful analysis of the NMR data may 
be needed: e.g. the power-law temperature dependence 
should be determined precisely to obtain the correct coef- 
ficient as a function of the magnetic field. Although such 
an enhancement has not beert-abserved yet for a typical 
ID spin-1/2 system CuHpCl,Ell we hope that more de- 
tailed experiments near the critical magnetic field may 
reveal the enhancement effect discussed here. 

Although we have studied a rather specific integrable 
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FIG. 8. The enhancement factor T± for x± around the 
cusp in Fig. [j](b). 
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model, some characteristic features found in this paper 
are expected to hold for more generic ID multicompo- 
nent quantum spin systems. In this connection, we note 
here that the SU(2) x SU(2) deformation of the mass- 
less SU(4) model drives the system to a massive phase 
with spin gaps for all thei-excitation modes, according 
to recent intensive studies. E3~EJ We think that the char- 
acteristic features obtained in the present study can be 
also applied to the static and dynamical properties in this 
model. 
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